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Abstract 

With an aim of establishing how the [S/Fe] ratios behave at the very low metallicity regime down 
to [Fe/H] ^ —3, we eonducted a non-LTE analysis of near-IR S i triplet lines (multiplet 3) at 10455- 
10459 A for a dozen of very metal-poor stars (—3.2 ^ [Fe/H] ^ —1.9) based on the new observational data 
obtained with IRCS-I-A0188 of the Subaru Telescope. It turned out that the resulting [S/Fe] values are 
only moderately supersolar at [S/Fe] ~ +0.2-0.5 irrespective of the metallicity. While this "flat" tendency 
is consistent with the trend recently corroborated by Spite et al. (2011, A&A, 528, A9) based on the 
S I 9212/9228/9237 lines (multiplet 1), it disaffirms the possibility of conspicuously large [S/Fe] (up to 
^ +0.8) at [Fe/H] ~ ~3 that we once suggested in our first report on the S abundances of disk/halo stars 
using S I 10455-10459 lines (Takeda & Takada-Hidai 2011, PASJ, 63, S537). Given these new observational 
facts, we withdraw our previous argument, since we consider that [S/Fe]'s of some most metal-poor objects 
were overestimated in that paper; the likely cause for this failure is also discussed. 
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1. Introduction 

The galactic evolution of sulfur (one of the a-group el- 
ements) has been a matter of controversy, since differ- 
ent [S/Fe] behaviors with a decrease of [Fe/H] were sug- 
gested in the regime of metal-poor halo stars depending 
on the lines used; i.e., ever- increasing [S/Fe] even up to 
~ +0.8 (S I 8693-4 lines of multiplet 6) or nearly constant 
[S/Fe] at a mildly supersolar value around ^ +0.3 (S i 
9212/9228/9237 lines of multiplet 1). See, e.g., Takeda et 
al. (2005) and the references therein for more details. 

Given this situation, we recently carried out a system- 
atic study on the [S/Fe] ratios of 33 disk/halo stars over 
a wide range of metallicity (—3.7 < [Fe/H] ^ +0.3) while 
newly exploiting the S i triplet lines at 10455-10459 A 
(multiplet 3) based on the near-IR spectra obtained with 
Subaru IRCS+A0188, which had barely been used be- 
fore (Takeda & Takada-Hidai 2011; hereinafter referred 
to as Paper I). Rather unexpectedly, while the the local 
plateau of [S/Fe] ~ +0.2-0.4 (flat trend) was confirmed at 
—2.5 ^ [Fe/H] ^ —1.5 in consistent with the tendency al- 
ready established from the S i 9212/9228/9237 lines (e.g., 
Nissen et al. 2007), we found a considerably large [S/Fe] 
ratio amounting to ~ +0.7-0.8 dex at very low metallicity 
([Fe/H] ^ —3), which apparently makes a puzzling discon- 



tinuity in the narrow interval of —3 ^ [Fe/H] ^ —2.5. If 
this trend is real, the chemical evolution of sulfur would 
have to be considered differently from other a elements 
generally showing a plateau at a mildly supersolar [a/Fe] 
over the halo metallicity range. 

Soon after we published Paper I, however. Spite et al. 
(2011) reported new results of their extensive study on the 
[S/Fe] of extremely metal-poor stars, which are markedly 
against our conclusion (i.e., considerably high [S/Fe] at 

[Fe/H] 3). That is, using the S i 9212/9228/9237 

lines based on the VLT/UVES data, they showed that 
[S/Fe] ratio is almost constant at mildly supersolar val- 
ues of ^ 0.2-0.5 over the very metal-poor regime of 
—3.5 ^ [Fe/H] < —2.5.^ How should we interpret this 
discordance? Do multiplet 1 lines (~ 0.92 yum) and mul- 
tiplet 3 lines (^ 1.05 /im) yield different S abundances at 
the extremely low-metallicity regime? 

Yet, we need to ascertain in the first place that the 
trend of high [S/Fe] at [Fe/H] ^ —3 we obtained in Paper 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



Jonsson et al. (2011) also recently investigated the sulfur abun- 
dances of halo giants of —2.5 ^ [Fe/H] ^ —1.5 with the S I 
10455-10459 triplet (which we used in Paper I) and the forbid- 
den [S l] line at 10821 A based on the VLT/CRIRES data. They 
then concluded that [S/Fe] favors a fiat trend at ~ -1-0.4 and a 
high [S/Fe] (i.e., rise or scatter at low [Fe/H]) is rather unlikely. 
However, since the [Fe/H] range they studied is limited to be- 
tween ~ —2.5 and ~ —1.5 (where we also derived a near-plateau 
[S/Fc]), their finding does not contradict the conclusion of Paper 
I. 
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I from the S I 10455-10459 lines universally exists for 
extremely metal-poor stars in general, since that conclu- 
sion was extracted from only a few objects based on the 
spectra of not-so-sufiicient quality. We thus decided to 
reinvestigate the [S/Fe] behavior for a larger sample (ca. 
a dozen objects) specifically confined to very metal poor 
stars (—3.2 < [Fe/H] ^ —1.9) based on our new observa- 
tions lately conducted again with Subaru IRCS-I-A0188. 
The purpose of this article is to report the outcome of this 
new analysis. 

2. Observational Data 

The near-IR spectroscopic observations were conducted 
on 2011 August 17 and 18 (UT) by using IRCS-hA0188 of 
the Subaru Telescope for the selected 13 very metal-poor 
stars, among which 4 are turn-off dwarfs and 9 are evolved 
giants. In addition, we also observed Vesta in order to get 
the sun-light reference spectrum with the same equipment. 
The list of the targets is given in table 1 . Note that G 64- 
37 (which was studied in Paper I) was again included in 
the present sample. The details of the instrument and 
its setting along with the data reduction procedures (the 
same as in our previous observations in 2009 July) are 
described in section 2 of Paper I. 

In this observing run, a special attention was paid to 
achieve a sufficiently high signal-to-noise ratio so as to de- 
tect very weak S i 10455-10459 lines of extremely metal- 
poor stars. For this purpose, especially long integrated 
exposure times (~ 1-4 hours) were expended to the com- 
paratively faint ( J ~ 9-11 mag) lowest metallicity stars 
([Fc/H] - -3), such as CS 30323-048, HD 126587, G 206- 
34, G 64-37, HE 1523-0901, and BD-16 251. For the 
other brighter objects of J ^ 5-8 mag, we set the expo- 
sure times from a few minutes to ^ 20 min depending on 
the brightness. The S/N ratios, eventually accomplished 
in the neighborhood of the S i lines in our finally resulting 
z J-band spectra (with the resolving power of i? ^ 20000), 
are typically - 300-500 (cf. table 1).^ 

3. Analysis and Results 

The atmospheric parameters {Tos, fogff, wt, and [Fe/H]) 
necessary for constructing the model atmosphere for each 
star were taken from various published studies (cf. table 
1). Then, as in Paper I, the S abundance was evaluated 
by way of the non-LTE spectrum-synthesis analysis while 
applying Takeda's (1995) automatic fitting procedure to 
the region of S i 10455-10459 lines. 

Since we modeled the observed stellar line profile (Dobs) 
by the convolution of (i) the intrinsic spectrum {Dq] where 
only the elemental abundance A is allowed to vary since 
the model atmosphere and the microturbulence are given) 



Actually, we noticed that some ripple patterns appeared in sev- 
eral restricted portion of the spectrum, which may be attributed 
to an imperfect flat-fielding. While this pattern was found to 
fall on the region of the S I lines in several cases depending on 
the stellar radial velocity, we could successfully remove them by 
dividing the spectrum by that of a rapid rotator. 
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Fig. 1. Synthetic spectrum fitting for the S I 10455-10459 
triplet lines. The best-fit theoretical spectra are shown by 
(blue) solid lines, while the observed data are plotted by (red) 
symbols. In each of the panels, spectra are arranged (from top 
to bottom, from left to right) in the ascending order of [Fe/H] 
as in table 1. A vertical offset of 0.05 is applied to each 
spectrum relative to the adjacent one, and the wavelength 
scale has been adjusted to the laboratory frame. 

and (ii) the Gaussian macro-broadening function /m(i') 
[oc exp(— ?;^/?;^)] parametrized by vm (including the com- 
bined effects of the instrumental broadening, the macro- 
turbulence, and the projected rotational velocity), such 
as I?obs = Dq * /m, adjustable free parameters in accom- 
plishing the best fit are A and vm , both of which were ac- 
tually varied in the previous analysis of Paper I. However, 
for the reason mentioned in the next section, we inten- 
tionally fixed the vu parameter in this S i 10455-10459 
fitting for five extremely metal-poor stars (CS 30323- 
048, HD 126587, G 206-34, G 64-37, HE 1523-0901, and 
BD— 16 251) at the pre-determined values which had been 
established in advance from the analysis applied to the 
strong C-I-Si feature at 1.069 iim. [where A(C), A{Si), 
and Vm were varied to search for the best fit]. The fi- 
nal Vm values resulting from (or assumed as fixed in) the 
profile- fitting analysis are given in table 1. 

How the theoretical spectrum for the converged solu- 
tions fits well with the observed spectrum is displayed in 
figure 1, and the resulting non-LTE S abundances (A^) 
along with the [S/Fe] values'^ are presented in table 1. 

In figure 2a are plotted the resulting [S/Fe] ratios 
against [Fe/H], where the abundance uncertainties (Stqv', 
cf. subsection 4.2 in Paper I) caused by ambiguities in 
Toff (±100 K), log5 (±0.2 dex), and Vt (±0.3 km s~^) 
are shown by thin error bars (though they are typi- 
cally on the order of ±0.1 dex and not very signifi- 
cant).'^ We also derived EW10455 (equivalent width in- 



[S/Fc] = [S/H] - [Fe/H], where [S/H] ee - 7.20. We used 

^0 " 



j45 = 7.20 (the value derived in Paper I based on the solar flux 



spectrum atlas of Kurucz et al. 1984) as the reference solar 
abundance in order to keep consistency with our previous study, 
which anyhow matches the present result (7.21; cf. table 1) 
derived from the spectrum of Vesta very well. 
We should remark that largely different atmospheric parame- 
ters from those adopted here are reported in the literature, es- 
pecially for high-gravity turn-off stars. For example, Ishigaki, 
Chiba, and Aoki (2010) derived (T<,fi, logg, [Fe/H]) of (5821, 
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Fig. 2. Sulfur abundances and the related quantities plotted 
against [Fe/H]: (a) [S/Fe] corresponding to non-LTE sulfur 
abundance, where attached thin error bars represent the am- 
biguities due to uncertainties in the atmospheric parameters 
i^Tgv'} rf. subsection 4.2 in Paper I), (b) A10455 (non-LTE 
correction for the S I 10455 line), (c) EW10455 (equivalent 
width for the S I 10455 line) with the S/N-dependent intrin- 
sic random error {±5EW) estimated from Cayrel's (1988) for- 
mula, (d) ^^^W'ci'r (abundance variation in response to EW 
changes of ±5EW). Dwarfs (logg > 3) and giants (logg < 3) 
are indicated by filled (scarlet) and open (brown) squares, 
respectively. 

vcrscly computed from A^) along with the corresponding 
non-LTE correction (A5^q455; negative values with typi- 
cal extents of ^ 0.1-0.3 dex) for the strongest component 
(at 10455.45 A) of the triplet in the same manner as de- 
scribed in subsection 4.2 of Paper I, which are given in 
table 1 and shown in figures 2b and 2c, respectively. The 
error bars (a few mA) accompanied with EWioi55 in fig- 



3.50, -3.66) and (6105, 3.65, -3.36) for CS 30323-048 and G 64- 
37, respectively, which are considerably discrepant from Nissen 
et al.'s (2007) results of (6338, 4.32, -3.21) and (6432, 4.24, 
—3.08) adopted in this study; i.e., by ~ 300-500 K lower in 
Tgff, by ~ 0.6—0.8 dex lower in logg, and by ~ 0.3—0.5 dex in 
[Fe/H]. Interestingly, since the signs of the abundance correc- 
tions are opposite (e.g., St— = -1-0.03 dex for ATcfi = —100 K, 
Sg- = -0.06 dex for Alogg = -0.2 dex, for G 64-37; cf. elec- 
tronic table E2 in Paper I), the net effect on A{S) itself is only of 
~ 0.1 dex level and thus not so significant even in such extreme 
cases. This means, however, that [S/Fe] would be appreciably 
raised because of the lowered [Fe/H], which eventually makes 
a serious discordance in comparison with the [S/Fe] results of 
other giants stars. 



ure 2c are the S/N-dependcnt uncertainties (SEW) esti- 
mated by Cayrel's (1988) formula, ^ 1.6{w SxY^^e, where 
w is the typical line FWHIVI 0.5-1 A; assumed to be 

0. 75 A), Sx is the pixel size (0.25 A), and e is the photon- 
statistics accuracy (^ (S/N)^^). The abundance errors 
^^EWerr^ in rcsponsc to these uncertainties in EW10455 
are further displayed in figure 2d, where we can see that 
these errors become considerable (e.g., a few tenths dex) 
for extremely metal-deficient stars ([Fe/H] ^ —3) showing 
very weak hues {EW10455 ^5-6 mA). 

4. Discussion 

We are now ready to answer the question which moti- 
vated this study: "How do the [S/Fe] ratios of very metal- 
poor stars behave around [Fe/H] ^ —3: a sudden jump to 
a considerably high value of ^ -1-0.8 (as suggested in Paper 
I) or only a mildly supersolar value at ^ -t-0.4 maintaining 
a flat trend (as derived by Spite ct al. 2011)?" We can 
recognize in figure 2a that the [S/Fe] ratios do not show 
any appreciable increase with a decrease in the metallic- 
ity, which are almost constant over —3.2 ^ [Fe/H] ^ —1.9 
with the mean of ([S/Fe]) ~ -1-0.34 (ct = 0.16) irrespective 
of dwarfs and giants. These [S/Fe] values of very metal- 
deficient stars resulting from this study are combined with 
those of halo/disk stars derived in Paper I on the [S/Fe] vs. 
[Fe/H] plot displayed in figure 3, from which we can state 
that [S/Fc] ratios arc almost "flat" around ^ +0.3 over 
the wide metallicity range between [Fe/H] ^ — 1 and ~' —3. 
Thus, we have to admit that Spite et al.'s (2011) argument 
actually represented the truth, which (in turn) means that 
our previous conclusion of Paper I was not correct, as far 
as the run of [S/Fe] around the regime of [Fe/H] ^ — 3 
is concerned. From an impartial point of view, however, 
this is a gratifying consequence in a sense, because the 
different abundance indicators (S i 9212/9228/9237 and 
S I 10455-10459) eventually turned out to yield consis- 
tent results with each other. 

Then, how should we interpret the prominently high 
[S /Fe] results once obtained for all three most metal-poor 
stars around [Fe/H] —3 among the sample in Paper I; 

1. e., ([S/Fe], [Fc/H]) = {+0.69, -3.08) [G 64-37], (-fO.76, 
-3.06) [BD-18 5550], and (-1-0.66, -2.71) [HD 115444]? 
With an intention to confirm whether or not these results 
are really reliable, we reexamined the process of the pre- 
vious analysis for these stars. 

Regarding BD-18 5550 and HD 115444, we noticed 
that the solutions of vm (macrobroadening parameter; cf. 
section 3) resulting as by-products of spectrum fitting arc 
unusually high (18.6 and 19.4 km s~^, respectively) com- 
pared to other stars where values around 10 km s~^ 
are in common.^ However, additional spectrum fitting 



We can roughly make an order-of-magnitude estimate of 
including three kinds of macrobroadening [instrumental profile 
broadening (?^ip), macroturbulence broadening {vmt), and ro- 
tational broadening (frt)] based on a simple modeling, where 
all three broadening functions are assumed to be Gaussian 
and each of these parameters are the corresponding e-folding 
width, between which the following relation holds: = 
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Fig. 3. [S/Fe] vs. [Fe/H] relation based on the results in ob- 
tained this study (larger squares) combined with those derived 
in Paper I (smaller circles). As in figure 2, open and filled 
symbols are for giants (logg < 3) and for dwarfs (logg > 3), 
respectively. 
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analyses to the C+Si feature at A 1.069 /im carried 
out for these stars yielded quite reasonable vyi values of 
9.4 km s-i (BD-18 5550) and 10.2 km s'^ (HD 115444) 
(of. figures 4a and 4c). Since almost the same wm should 
(in principle) result from different lines, we can not help 
considering that the wm values for these two stars obtained 
from the S I 10455-10459 fitting in Paper I were inade- 
quately overestimated, which we suspect may presumably 
related to the profile-fitting technique we adopted. That 
is, the automatic solution-search algorithm (Takeda 1995), 
which simultaneously varies several parameters to accom- 
plish the best fit, is very efficient in case where the stellar 
line profiles are well defined. However, we should be care- 
ful when it is applied to the very weak-line case where 
noises are comparable to the signal of stellar lines, since it 
may yield physically meaningless solutions where profiles 
are appreciably damaged by noises. In such cases, vyi had 
better be fixed at a (more reliable) value derived from 
other stronger lines (such as the C-|-Si feature), instead 
of varying both A(S) and I'm, which is the approach we 
adopted for the five most metal-poor stars in this study 
(cf. section 3).^ Accordingly, we redetermined the S 

v^^ -f iJ^j. -f v^^ . Under this approximation, tijp is evaluated as 

^ip = FWHM/(2Vln2) ~ 9 km s"! (FWHM ~ 15 km s"! for 
R~ 20000). Then, we may expect that iimt would be as small 
as ~ 2 km s~^ according to the relation ^!mt — 0.4fj^x (cf. foot- 
note 12 of Takeda et al. 2008) since the typical fj^x (radial- 
tangential macroturbulence) is ~ 5 km s"'^ for early G dwarfs 
and early-K giants (cf. figure 17.10 in Gray 2005). Finally, we 
may regard that v-^t (— 0.94i)osini; cf. footnote 12 of Takeda 
et al. 2008) is of minor importance (e.g., presumably no larger 
than ~ 2—3 km s^^ in most cases), since the stellar rotation 
must have been spun down in these old halo stars. Accordingly, 
vyi is reasonably expected to be around ~ 10 km s^^, because 
it is primarily determined by the contribution from -Uip. 
® For the other stars, both A(S) and dm were varied to determine 
as in Paper I. However, we cheeked for each star that the result- 
ing «m(S) (cf. table 1) is not in serious disagreement with -um 
(C-|-Si). Note that HD 195636 is a rather unusual star, since 



Fig. 4. Reanalysis of two stars (BD-18 5550 and HD 115444 
in the left and right panels, respectively) which showed partic- 
ularly high [S/Fe] values in Paper I, in order to demonstrate 
how the choice of vyy (macrobroadening parameter) influences 
the resulting sulfur abundance derived from spectrum fitting. 
Upper panels (a) and (c) ■ ■ ■ Derivation of best-fit dm from 
the C+Si feature at ~ 1.069 /im. Lower panels (b) and (d) 
• ■ ■ Determination of A(S) with two kinds of dm treatments, 
where dm is varied as an adjustable parameter (upper spec- 
trum; method adopted in Paper I) or dm is fixed at the value 
derived from the C-|-Si feature (lower spectrum). 

abundances from the S i 10455-10459 lines while fixing 
the wm values at the values derived from the C+Si fea- 
ture, and obtained A(S) values (4.67 and 4.83), which 
are appreciably lower than the results in Paper I (4.90 
and 5.15) by 0.23 dex and 0.32 dex for BD-18 5550 and 
HD 115444, respectively (cf. figures 4b and 4d). We thus 
consider that the [S/Fe] values of these two stars reported 
in Paper I should be revised downward by these amounts, 
which eventually makes -hO.53 (BD-18 5550) and -hO.34 
(HD 115444); i.e., being almost within the [S/Fe] range 
concluded in this study. 

The star G 64-37 was actually an object of special 
attention, since the S abundance (based on S i 10455- 
10459) derived in Paper I turned out appreciably higher 
(by ~ 0.6 dex) than Nissen et al.'s (2007) result (based on 
S I 9212/9228/9237) only for this extremely metal-poor 
star, despite that a good agreement is seen for other com- 
paratively more metal-rich stars; this fact motivated us to 
include this star again in our target list. Unlike the case 
of BD-18 5550 and HD 115444, a serious mismatch of 
wm was not found in the analysis of Paper I. Instead, we 
consider that the essential problem was the comparatively 
poor S/N ratio 100) of the spectrum used in Paper I, 
which is evidently insufficient for a reliable abundance de- 
its vyi (17.0 km s~^) is appreciably large (an alternative C-|-Si 
fitting also yielded a similar result), which might have a com- 
paratively high Desinj for stars of this class. 
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Fig. 5. Comparison of the spectra and the A(S) determi- 
nation procedures for G 64-37, for which two independent 
analyses were performed in Paper I (based on the 2009 data) 
as well as in this study (based on the 2011 data). 

termination from very weak lines with depression on the 
order of ~ 1%. The observed and the fitted spectrum in 
Paper I and those in this study are compared in figure 
5. We can see from this figure that the previous analy- 
sis was severely influenced by the large noise, while the 
situation has been improved in present case (S/N ^ 200), 
resulting in an appreciable difference of solution between 
the two cases. We naturally consider the present result of 
A{S)= 4.54 ([S/Fe] = +0A2) is more credible than that 
in Paper I, requiring a downward revision of the previ- 
ous [S/Fe] value (+0.69) by -0.27 dex, by which the dis- 
agreement with Nissen et al. (2007) has been reasonably 
mitigated by this amount. 

5. Conclusion 

This study was motivated by the recent extensive work 
done by Spite et al. (20ii) who reported that [S/Fe] 
ratios of very metal-poor stars determined from the S i 
9212/9228/9237 lines show a nearly flat behavior at a 
mildly supersolar level of ~ +0.2-0.5 over a wide metal- 
licity range of —3.5 ^ [Fe/H] ^ — 2, which markedly con- 
tradict the conclusion of our previous study based on S i 



10455-10459 lines (Paper I) that such a flat tendency of 
[S/Fe] (persisting down to [Fe/H] —2.5) is followed by 
a sudden jump up to -1-0.7-0.8 around [Fe/H] ~ —3. 

With an intention to resolve the cause of this discrep- 
ancy, we rechallenged the task of clarifying the [S/Fc] 
ratios at the extremely low metallicity regime down to 
[Fe/H] ~ —3 by using the same triplet lines as used in 
Paper I, based on the new observational data for an ex- 
tended sample of 13 very metal-poor stars observed with 
IRCS-hA0188 of the Subaru Telescope. 

In almost the same manner as in Paper I, we con- 
ducted a non-LTE spectrum fitting analysis of S i 10455- 
10459 triplet, and found that the resulting [S/Fc] values 
were moderately supersolar uniformly scattering around 
- +0.3-0.4 [with the mean abundance of ([S/Fe]} = +0.34 
(cr = 0.16)] over — 3.2 < [Fe/H] < -1.9 without any sys- 
tematic [Fc/H]-depcndencc, which confirmed the conse- 
quence corroborated by Spite et al. (2011) based on the 
S I 9212/9228/9237 lines. 

Given these new observational facts, we reexamined 
the process of our previous analysis for the three ex- 
tremely metal-poor stars (G 64-37, BD— 18 5550, and 
HD 115444), for which we derived prominently high [S/Fe] 
values +0.7-0.8) that eventually lead to the conclu- 
sion of Paper I. Regarding G64-37, our reanalysis using a 
new spectrum of higher quality yielded a result lower than 
the previous value by ~ 0.3 dex. For BD— 18 5550 and 
HD 115444, we noticed that our automatic profile-fitting 
method (which varies both the abundance and the broad- 
ening parameter to find the best fit) resulted in unreason- 
ably large solutions of the broadening width, because of 
the considerable weakness of the line profile severely dam- 
aged by noises. When the broadening parameter was fixed 
at the more reasonable values determined from stronger 
lines, we found that the revised solution of the S abun- 
dance is lowered by ^ 0.3-0.4 dex for both of these stars. 
Consequently, it is likely that we had overestimated the 
[S/Fe] values of these three stars in Paper I by ^ 0.3- 
0.4 dex. 

Accordingly, we now consider that the flat trend of 
[S/Fe] (without any systematic rise with a decrease of 
metallicity) represents the truth, at least with regard to 
the overall [S/Fe] behavior of very metal-poor stars in 
general, which means the withdrawal of our previous argu- 
ment that [S/Fe] experiences a sudden jump up to conspic- 
uously large [S/Fe] of ^ +0.8 as the metallicity is lowered 
down to the [Fe/H] ~ — 3 regime. 

We flnally remark, however, that the existence of some 
stars deviating from the main trend with appreciably 
higher/lower [S/Fe] amounting to [S/Fc] ~ +0.7/0.0 as 
a result of the natural diversity is not necessarily be ex- 
cluded (e.g., we derived [S/Fe] = +0.66 for HD 126587 and 
[S/Fe] = -0.01 for HD 13979). Actually, since our 13 stars 
(at [Fe/H] < -2) yielded ([S/Fe]) = +0.34 with a (stan- 
dard deviation) of 0.16, we may expect the probability of 
finding > la deviation ([S/Fe] < +0.2 or [S/Fe] > +0.5) 
and > 2(7 deviation ([S/Fe] < 0.0 or [S/Fe] > +0.7) to be 
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~ 30% and ^ 5%, respectively/ Form this point of view, 
the recent Koch and Caffau's (2011) result of [S/Fe] as 
somewhat high as ^ +0.5 for a red giant in the very metal- 
poor globular cluster NGC 6397 ([Fe/H] 2.1) may be 

understandable without invoking a bimodal [S/Fe] trend 
such as they discussed. 
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This estimation naturally depends on the extent of cr. We note 
that the fT([S/Fe]) of [Fe/H] < —2 stars derived by other inves- 
tigators are somewhat smaller than our result; e.g., cr = 0.07 
(Nissen et al. 2007), rr = 0.12 (Spite et al. 2011), cr = 0.11 
(Jonsson et al. 2011). It should thus be important to quantita- 
tively establish not only the trend of [S/Fe] on the average but 
also the extent of its dispersion based on a large sample of very 
metal-poor stars. 
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Table 1. Parameters of the program stars and the results of abundance analyses. 



Name 




logs 


vt 


[Fe/H] 


Ref. 


vm 


A 


-EVV^10455 


Al0455 


[S/Fe] 


S/N 


Remark 






(K) 


(cm s~^) 


(km s"^) 


(dcx) 




(km s"^) 


(dox) 


(mA) 


(dox) 


(dex) 








CS 30323-048 


6338 


4.32 


1.5 


— 3.21 


NIS07 


(8.7) 


4.29 


2.8 


— 0.25 


-1-0.30 


350 


I'M fixed, larf 


yer uncertainty 


HD 126587 


4700 


1.05 


1.7 


—3.16 


HAN 11 


(12.7) 


4.70 


14.2 


— 0.27 


4-0.66 


300 


VM fixed 




G 206-34 


5825 


3.99 


1.5 


-3.12 


RIC09 


(9.4) 


4.55 


3.6 


-0.19 


H-0.47 


450 


I'M fixed 




G 64-37 


6432 


4.24 


1.5 


-3.08 


NIS07 


(8.3) 


4.54 


5.5 


-0.26 


-1-0.42 


200 


vm fixed 




HE 1523-0901 


4630 


1.00 


2.6 


-2.95 


FRE07 


(10.3) 


4.42 


6.7 


-0.19 


4-0.17 


350 


vm fixed 




BD-16 251 


4825 


1.50 


1.8 


-2.91 


ANDIO 


9.4 


4.64 


10.2 


-0.23 


-1-0.35 


450 






HD 195636 


5370 


2.40 


1.5 


-2.77 


CAR03 


17.0 


4.88 


19.1 


-0.35 


-1-0.45 


300 


large i?M (cf. 


footnote 5) 


G 186-26 


6417 


4.42 


1.5 


-2.54 


NIS07 


12.0 


5.02 


10.7 


-0.15 


+0.36 


500 






HD 186478 


4730 


1.50 


1.8 


-2.42 


HANll 


13.3 


5.19 


22.8 


-0.19 


4-0.41 


350 






HD 6268 


4735 


1.61 


2.1 


-2.30 


SAI09 


9.5 


5.20 


21.4 


-0.17 


4-0.30 


450 






HD 13979 


5075 


1.90 


1.3 


-2.26 


BUROO 


11.8 


4.93 


16.6 


-0.21 


-0.01 


450 






HD 221170 


4560 


1.37 


1.6 


-2.00 


SAI09 


9.2 


5.49 


28.7 


-0.16 


4-0.29 


200 






HD 216143 


4525 


1.77 


1.9 


-1.92 


SAI09 


9.4 


5.50 


19.9 


-0.10 


4-0.22 


550 






Vesta (Sun) 


5780 


4.44 


1.0 


0.00 




10.7 


7.21 


123.5 


-0.09 


4-0.01 


600 







In columns 1 through 6 are given the star designation, effective temperature, logarithmic surface gravity, microturbulent velocity dispersion, Fe 
abundance relative to the Sun, and key for the reference of atmospheric parameters: BUROO ■•• Burris et al. (2000), CAR03 ■■■ Carney et al. 
(2003), FRE07 ■ ■ ■ Frobel et al. (2007), HANll ■ ■ ■ Hansen and Primas (2011), NIS07 ■ ■ ■ Nissen et al. (2007), RIC09 ■ ■ ■ Rich and Bocsgaard (2009), 
SAI09 ■ ■ ■ Saito et al. (2009). Columns 7-11 present the results of the abundance analysis based on the S I 10455-10459 profile-fit: vm is the best-fit 
macrobroadening parameter (while those in parentheses are the assumed or fixed values, which were separately derived from the C-^Si 1.069 fim 
feature fitting), is the non-LTE logarithmic abundance of S (in the usual normalization of H — 12.00) derived from spectrum-synthesis fitting, 
-EW^i0455 is the equivalent width (in mA) for the S I 10455 line inversely computed from , A10455 is the non-LTE correction (= — ^10455) 
for the S I 10455 Une, and [S/Fe] (= — 7.20 — [Fe/H]) is the S-to-Fe logarithmic abundance ratio relative to the Sun. The S/N ratio of the 
spectrum estimated at the position of the S I lines is given in column 12. The objects are arranged in the ascending order of [Fe/H]. 



